One-dimensional Surface Bound States in d-wave Superconductors 
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We present an exact quantum theory for the bound states in the vicinity of an edge or a line of 
impurities in a d x 2_ y 2 superconductor. For a (llO)-surface we show that a finite dispersion of the 
one dimensional band of bound states leads to a two peak structure in the density of states (DOS). 
We study the effect of an applied magnetic field and a subdominant id xy order parameter on the 
DOS and discuss the implications of our results for tunneling experiments. 



PACS numbers: 73.20.-r, 73.20.At, 74.25.Jb 

The existence of a zero bias conductance peak 
(ZBCP), observed by in-plane tunneling spectroscopy 
on YBa 2 Cu 3 7 (YBCO) thin films @-|| has been in- 
terpreted as a signature of unconventional superconduc- 
tivity. In the framework of the quasiclassical approxi- 
mation, it was first shown by Hu [p|, and subsequently 
studied in more details in Refs. (^,0), that in supercon- 
ductors with d x 2_ y 2 symmetry the ZBCP arises from 
a one-dimensional dispersionless band of surface bound 
states. Recently Covington et al. discovered the split- 
ting of the ZBCP into two peaks not only in the pres- 
ence of an applied field, H, but also for H = below 
T s = 8 K |. Within the quasiclassical approach, the 
splitting for H = was ascribed to a spontaneously bro- 
ken time reversal symmetry (BTRS) on the surface of 
YBCO and the generation of a subleading superconduct- 
ing order parameter with id xy or is symmetry Q. While 
in conventional superconductors corrections to the bound 
state energy beyond the quasiclassical approximation are 
too small to be observed experimentally, the corrections 
in high temperature superconductors (HTSC) are of the 
order of A^/Ep 3 meV and thus comparable to the ex- 
perimentally observed peak splitting. This necessitates a 
fully microscopic treatment of the problem before a con- 
clusion about the nature of the splitting for H = can 
be reached. 

In this Letter we present a microscopic theory which 
goes beyond the quasiclassical approximation and allows 
us to capture properties of the ZBCP that eluded earlier 
investigations. By developing a scattering matrix for- 
malism for one-dimensional (ID) defects, appropriate to 
describe surfaces or lines of impurities, we show that the 
band of bound states has a finite dispersion that leads to 
a two-peak structure in the density of states (DOS), even 
in the absence of a BTRS. We predict the existence of 
Friedel oscillations associated with the bound states and 
show that their spatial dependence varies strongly with 
the so-called acceptance angle of the measurements. In 
the presence of a supercurrent or magnetic field the DOS 
exhibits a four-peak structure for small acceptance an- 
gles, but only two peaks with increased frequency split- 
ting for large acceptance angles. We show that the onset 
of a subdominant order parameter with id xy symmetry 
leads to a sharpening of the two-peak structure in the 



DOS, but does not induce any qualitatively new features. 
Finally, we compare our results with the in-plane tunnel- 
ing experiments on YBCO || and comment on the ex- 
tensions of our formalism to treat surface bound states 
in systems with other non-trivial order parameters. 

We begin by presenting an exact scattering theory 
which describes the emergence of a bound state in the 
vicinity of an interface, e.g., a line of impurities or a 
surface along the (ImO) direction, in a d x 2_ y 2 supercon- 
ductor. The scattering of the electrons off the interface 
is described by U(r) = U J2r H r ~ Ri)> wnere the sum 
runs over all R; = (x, mx). The Green's function in the 
presence of this potential is given by 

G(r,r',w) = G (r-r,w) + 

; G (r - s, u)t (s, s', w)G (s' - r', w) , (1) 



where the fourier transform of the scattering T-matrix 
only depends on the momentum parallel to the line, q\\, 



f(q l{) uj) = [l-Uf 3 



' dp 

2^ 



(2) 



T a are the Pauli matrices, and the prime restricts the 
integration to those momenta in the first Brillouin zone 
(BZ) for which p x + mp y — q\\ = 2wn, with n integer. For 
the Green's function of the clean system one has 



(3) 



G (k,w) = [u,f - e k f 3 - A k 0) n]" 1 , 

where A k °^ = Ao(cos/c2; — cosfcj,)/2 is the d x 2_ y 2 super- 
conducting gap . For the normal state we use a quasi- 
particle dispersion, representative for the HTSC 

ek = — 2i(cos k x + cos k y ) — At' cos k x cos k y — /i (4) 

with t = 300 meV, t' /t = -0.3, fi/t = -0.99 and A = 25 
meV 0. 

In the following we consider for definiteness a surface 
along the (110) direction (with U — oo). The induced 
surface bound states represent quasiparticles that are lo- 
calized in the direction perpendicular to the surface, but 
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FIG. 1. Dispersion, oj s (q\\), (solid line) and gap edge of the 
particle-hole continuum (dashed line) of the surface bound 
states as a function of q\\ . Inset: characteristic FS of the 
HTSC and lines of integration for a (llO)-surface. We set the 
lattice constant ao = 1. 

have a well defined momentum parallel to it. The disper- 
sion, u s (qu), of these bound states is determined by the 
poles of the T-matrix and thus follows from the condition 

/ d * V-?-A2 =± J dp ^_" 2 S _ A2 ■ ( 5 ) 

J ^ s L p p J s L p p 

For a (llO)-surface, the integration is restricted to a line 
of momenta parallel to the (llO)-direction (dashed line 
in the inset of Fig. Q). For small q\\ the integrals in Eq.(||) 
can be expanded by standard means, and one finds that 
the dispersion is quadratic in q\\ if the Fermi surface (FS) 
is particle-hole symmetric and linear in qn when this sym- 
metry is absent (as is the case for the FS shown in Fig. [I]). 
In Fig. [I] we present the dispersion of the bound states, 
obtained from the full numerical solution of Eq.(|5|), to- 
gether with the gap edge of the particle-hole continuum. 
The sharp drop of the bound state energy at q?, occurs 

when the line of integration crosses the FS at k c (see in- 
set). Just below this crossing, i.e., for q\\ < q?,, the poles 

of the T-matrix lie within the particle-hole continuum, 
and the surface states are damped resonant state. 

The DOS, renormalized by the defect/surface scatter- 
ing, at a distance r± = (r x — r y )/\/2 from the surface 
is given by N(ui,r±) = Im [G (r, r, lu) + SG(r, r, oj)] u , 
where 

6d(r ± ,u)= [ A p-&i{u,q h T)t(q h u)& (u,q h -T) 
J-A 2vr 

G' { " ) (u,q h r)= J |E e ^G (^p), (6) 

and the double prime restricts the integration to all mo- 
menta in the first BZ Witli p x +mp y ~ q\\ — 0. Summation 
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FIG. 2. N(r±,u>) for (a) A = l/y/2, and (b) A = ^/2ty. 
The coloring indicates the magnitude of N(r±,tj). 

over all gy in Eq.(|^) corresponds to A = \/2tt. However, 
in the analysis of tunneling experiments it is in general 
assumed that A < \/2it due to the finite acceptance angle 
of the experiments ; in what follows we therefore con- 
sider two different values of A in our calculations of the 
DOS. Note that the above expression for 5G is similar to 
the one for a single impurity problem flT|| . 

In Fig. we plot N(r±,uj) obtained numerically for 
A = 1/V2 and \f2-n. In both cases, the DOS ex- 
hibits a two-peak structure and Friedel oscillations (FO) 
whose spatial dependence is described by N(r±,oj) ~ 
e lk ° r± - e~ r± ^. Both ko and £ vary with frequency and 
through ui s (qii) implicitly depend on qn. For A = 1/v2j 
there exist only two bound states with wave-vector ±gy 
for any frequency ui < to s (A) and one has ko(uj s ) = 
2kp.±(u) s ) = y/2(kp — kp) (see Fig. |l|). Moreover, since fco 
increases with increasing the wavelength of the FO 
decreases with increasing \uj\, as can clearly be seen in 
Fig. |^a. The lengthscale, £, for the exponential decay of 
the FO_into the bulk is set by £ = v F ^/^/A 2 - w| > 1, 
where A is the superconducting gap at the respective FS 
crossing. For A = \Z2tt, where we integrate over all gn, 
the FO below and above 4 meV are qualitatively different. 
For uj < 4 meV there exist four bound states, two close 
to <7|| = and q^ = it (see Fig. ^), whose superposition 
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FIG. 3. N(u),r±) for p s = 0.015 and (a) A = l/y/2, and 
(b) A = y/2-rr 

leads to the rapid decay of the FO at low frequencies. 
In contrast, for u> > 4 meV, only two bound states exist 
for a given frequency and the DOS again exhibits well 
defined FO. While a particle- hole symmetry is absent in 
the DOS for a given r±, it is recovered after averaging 
over r±. 

We next study the effects of an applied magnetic field 
on the DOS. To avoid the problems arising from the 
presence of vortices and a spatial variation of the su- 
percurrent, we consider for simplicity an average super- 
current with momentum p s ; our results are thus exact 
for the case when a uniform supercurrent is applied to 
the sample. For p s ^ the electronic spectrum experi- 
ences a Doppler-shift dD E k = ^Jt\ + |A k | 2 +Vi?(k)-p s , 
where Vi?(k) = de^/dk. Thus, by making the substitu- 
tion lo — > u — vp(k) ■ p s in Eq.([y), we obtain N(u, rj_) 
which we plot in Fig. || for a supercurrent along the (1 In- 
direction with p s = 0.015. The bound state energy is now 
shifted to ui s {<1\\) = w s(<Z||) +vf ■ Ps, where for ±q\\, vj? 
is the Fermi velocity at the respective FS crossing. This 
shift leads to a four-peak structure in the DOS, as can 
be seen in Fig. 3a for A = 1/ v2. Since the energy shift 
depends on qu through Vp, the peaks which are shifted 
towards u> = are sharpened, but still possess the same 
spectral weight as the apparently weaker peaks shifted 
to larger \oj\. In contrast, for A = \/2ir (Fig. ^d), the 
Doppler-shift varies much more strongly along the lines 
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FIG. 4. N(u, r±) as a function of u> for \/2r± = 1.25 for a a 
superconducting gap with d x 2_ y 2 (solid line) and d x 2_ y 2+id xy 
symmetry (dashed line). 

of integration which preserves the two-peak structure in 
the DOS, increases the frequency splitting between the 
two peaks, and leads to a filling-in of the gap between the 
peaks. Thus, we predict that only for small acceptance 
angles can a four-peak structure be resolved. 

We now turn to the effect of a subdominant or- 
der parameter with d xy symmetry and gap = 
Aisin/cz sinky on the DOS. In Fig. |we plot N(u>, r±) 
for a superconducting gap with d x 2_ y 2 (solid line) and 
d x 2_ y 2 + id xy symmetry (dashed line). For the latter, we 
assume that the total superconducting gap is given by 

A k 0) +* A k ) with A i/ A o = 0.1 g|. Thc D0S at sma11 
u) increases linearly with energy for the d x 2_ y 2 case and 
exhibits as expected a gap for the d x 2_ y 2 + id xy case. 
Note, that the the onset of an id xy component shifts the 
two peaks to only slightly higher frequencies. 

We next comment on the differences between our scat- 
tering matrix approach and earlier quasiclassical ap- 
proaches within the Bogolyubov-deGennes formalism. In 
the latter, one considers only quasiparticles with a well 
defined momentum on the Fermi surface, such that for a 
(llO)-surface the incident and reflected quasiparticles ex- 
perience a superconducting gap of equal magnitude and 
opposite sign, which leads to a zero energy bound states 
for all <7y This result is exact for gy = 0, since here 

A(k) = 0, but for any finite qu, the scattering between 
states away from the FS introduces corrections to the 
quasiclassical energy. Hence, within our scattering ma- 
trix formalism, the quasiclassical results are recovered by 
restricting the integration in Eq. (JsJ) to the Fermi surface 
which yields uj s = for all q\\ . The appearance of a 
dispersing band of bound states is thus a pure quantum 
effect |Ti|] ) which cannot be obtained by a quasiclassical 
approximation. 

Experiments of || on YBCO materials show a sin- 
gle peak above T s — 8 K <C T c , and the emergence of a 
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two-peak structure at lower temperatures. This and the 
small increase in the overall width of the ZBCP was at- 
tributed to the appearance of a subleading order param- 
eter (id xy or is) below T s . However, the observed broad- 
ening of the ZBCP below T s is much smaller than the 
splitting between the peaks. This is expected if the peak 
splitting arises not from the appearance of a BTRS, but 
from the finite bound state dispersion discussed above. 
A subleading order parameter may sharpen the two-peak 
structure at lower temperatures and suppress the low en- 
ergy part of the spectrum but is not predominantly re- 
sponsible for the peak splitting. Moreover, recent experi- 
ments on Bi2Sr2Ca208+,5 |l5| ], which find that the width 
of the ZBCP does not increase with the emergence of a 
two-peak structure in the DOS, are difficult to explain 
within the BTRS scenario. Currently, experiments using 
electron-spin resonance |]l6f and phase sensitive measure- 
ments Jl7t are under way to determine whether a BTRS 
exists on the surface of HTSC. In addition, detailed mea- 
surements of the low energy part of the DOS also provide 
important information since the DOS for a d x 2_ y 2 super- 
conductor is linear at low frequencies, whereas that for 
d + id or d + is order parameters is fully gapped. Another 
interesting feature of the experiments |8| is that only two 
peaks are observed in an applied magnetic field, which 
shift to higher energies with increasing field. One possi- 
ble explanation of this result is that the acceptance angle 
of the experiment (see Fig. |^b) is large, in which case 
the regime of smaller acceptance angles can be accessed 
in tunneling experiments by using stronger barriers be- 
tween the superconducting and normal materials. In this 
limit, we predict that an applied magnetic field induces a 
four-peak structure in the DOS (Fig. ||a). An alternative 
explanation is that the acceptance angle is small but that 
the magnetic field generates a large id (or is) component 
at the boundary that pushes the peaks to higher energies 
(Fig. U). Each of the original peaks should then be split 
into two, however such splitting may be too small to be 
observed experimentally. The latter scenario is consis- 
tent with a strong decrease in the DOS at zero energy 
with increasing magnetic field ^] . 

Finally, the bound states discussed above provide an 
interesting new class of one dimensional fermions whose 
dispersion can be changed continuously by applying a 
magnetic field (or supercurrent), with the possibility of 
making all fermions chiral |l8[ ]. Moreover, the formal- 
ism presented above can be generalized to study surface 
bound states in systems with other non-trivial order pa- 
rameters, such as d-density or spin density wave states 
suggested recently [[[9) . We find that a zero energy bound 
states also exists for a d-density wave state in the vicinity 
of a (110) surface. This bound state can be observed in 
tunneling experiments and may thus prove to be an im- 
portant feature of the underdoped cuprates; an account 
of our results will appear elsewhere [jl8| . 

In summary we present an exact scattering theory for 
bound states which are induced close to lines of impu- 
rities or surfaces in a d x i_ v ? superconductors. We find 



that these bound states disperse linearly with momentum 
<gf II , a dispersion that gives rise to a two-peak structure in 
the density of states. For small acceptance angles, we 
predict that an applied magnetic field leads to a split- 
ting of the two peaks into four, while for a large accep- 
tance angle, there exist only two peaks. Our results are 
consistent with recent in-plane tunneling experiments on 
YBCO and provide important insight into their interpre- 
tation. We present detailed predictions for the spatial 
dependence of the Friedel oscillations in the DOS that 
may be measured by the STM experiments near the sur- 
faces, grain boundaries, or lines of impurities. 
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